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基于 景观 格局 的 丘陵 山区 耕地 生态 风险 评价 
一 一 以 重庆 市 潼 南 区 为 例 
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摘 要 : [目的 /意义 ] 丘陵 山区 农田 宜 机 化 整治 会 改变 耕地 景观 格局 、 高 程 、 坡 度 、 微 地 貌 等 ,评价 其 生态 风险 
为 后 续 整 治 工作 提供 数据 参考 具有 重要 意义 。 本 研究 的 目的 为 评价 丘陵 山区 农田 宜 机 化 整治 对 耕地 生态 风险 的 改 
变 情 况 以 及 探究 生态 风险 与 耕地 坡度 之 间 的 关系 。[ 方 法] 以 重庆 市 济南 区 20 个 县 为 评价 单元 ， 基 于 2010 年 和 
2020 年 土地 利用 数据 ， 采 用 ArcGIS 10.8 和 Excel 软 件 计 算 景 观 格局 指数 ， 通 过 灶 权 法 确定 各 指数 的 权重 并 构建 生 
态 风险 评价 模型 ， 揭 示 生 态 风 险 时 序 空间 变化 特征 ; 基于 数理 统计 原理 ， 对 生态 风险 与 坡度 进行 相关 性 分 析 ， 探 


究 生 态 风险 与 坡度 的 关系 。[ 结 果 和 讨论 ] 2010 年 和 2020 年 两 个 时 期 ， 干 扰 度 指 数 由 0.97 下 降 为 0.94， 耕 地 整体 
抗 干扰 能 力 增强 ; 脆弱 度 指数 由 2.96 增 加 为 3.20， 耕 地 结构 更 加 脆弱 ; 生态 风险 值 由 3.10 下 降 为 3.01， 耕 地 生态 


态 风险 值 呈 极 显著 相关 。 农 田 宜 机 化 整治 应 重点 关注 潼 南 区 南部 


A 
安全 性 提高 。 两 个 时 期 生态 风险 区 域 主要 以 低 风险 区 和 较 低 风险 区 为 主 ， 低 风险 区 面积 增加 6.44% ， 较 低 风险 区 
面积 增加 6.17%， 中 风险 区 面积 增加 24.4%， 较 高 风险 区 面积 减少 60.70%， 高 风险 区 面积 增加 16.30%， 耕 地 生态 
安全 区 域 相对 增加 。 赫 地 坡度 主要 以 2 一 25 为 主 ， 耕 地 坡度 小 于 15 "时 坡度 面积 占 比 与 生态 风险 值 呈 负 相 关 ， 耕 
地 坡度 大 于 15° 时 坡度 面积 占 比 与 生态 风险 值 呈正 相关 关系 ,坡度 处 于 5°~8”、15 一 25 、25 以 上 时 坡度 面积 与 生 


区 域 ， 并 集中 于 耕地 坡度 处 于 5 ~8° 和 15 ~25° 


区 域 。[ 结 论 ] 通过 评价 济南 区 农田 宜 机 化 整治 前 后 


th 


也 生态 风险 并 分 析 生态 风险 与 耕地 坡度 的 相关 性 ， 表 明 农田 


化 整治 工作 应 重点 关注 耕地 坡度 处 于 5 ~8° 和 15 ~25° 区域。 
关键 词 : 农田 宜 机 化 整治 ; 耕地 ; 生态 风险 评价 ; 景观 格局 
中 图 分 类 号 : X826 文献 标志 码 : A 


宜 机 化 整治 可 以 降低 耕地 生态 风险 ， 耕 地 坡度 面积 占 比 可 作为 精准 指导 农田 宜 机 化 整治 的 重要 依据 ， 漳 南 区 宜 机 
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051 言 小 且 分 散 等 自然 条 件 制约 着 中 国 丘 陵 山区 农机 化 发 


A 


丘陵 山区 是 中 国 重要 的 粮油 、 果 蔬 和 茶叶 等 农 
产品 生产 供给 基地 。 然 而 ， 坡 地 多 、 坡 度 大 、 地 块 


展 ， 对 农业 经 济 发 展 造 成 限制 。 中 国 自 20 世 纪 90 
年 代 中 期 开展 丘陵 山区 土地 整治 ， 以 增加 有 效 耕 地 


面积 "。 经 过 多 年 的 探索 与 实践 ，2014 年 重庆 市 延 
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伸 土 地 整治 内 容 ， 开 展 农 田 宜 机 化 整治 ， 大 力 推广 
选择 道路 通达 且 坡 度 不 大 于 25* 的 耕地 ， 将 零碎 分 
散 的 小 地 块 改造 成 坡 式 梯田 、 缓 坡地、 水平 条 田 等 
大 地 块 了 。 然 而 ， 农 田 宜 机 化 整治 改变 了 耕地 景观 
原始 的 结构 ， 使 景观 分 割 、 破 碎 程 度 加 剧 ， 一 定 程 
度 上 影响 农田 生态 系统 的 稳定 性 o ERER TF, 
研究 耕地 生态 风险 可 以 为 农田 宜 机 化 整治 的 风险 防 
探 提供 决策 依据 ， 从 而 推进 丘陵 山区 农业 高 质量 
发 展 。 

目前 ， 生 态 风险 评价 主要 关注 土地 利用 变 
AG) 、 城 市 、 流 域 " 9、 土 地 整治 “等 方面 ， 
主要 评估 方法 有 相对 风险 模型 (Relative Risk Mod- 
el, RRM) 和 生态 风险 指数 法 (Ecological Risk In- 
dex, ERI) “2 。 相 对 风险 模型 用 于 多 个 生态 系统 
的 区 域 复合 压力 生态 风险 评价 ， 如 陆地 、 水 域 等 ; 
生态 风险 指数 法 通过 选择 与 生态 风险 值 相关 的 景观 
格局 指数 构建 评价 模型 ， 用 于 景观 格局 的 研究 。 现 
有 基于 景观 格局 的 土地 整治 生态 风险 评价 主要 为 评 
价 指 标的 选取 以 及 生态 风险 评价 模型 的 构建 。 张 正 
峰 等 ”选取 田 块 规模 、 田 块 分 维 数 、 田 块 分 离 度 
等 景观 格局 指数 评估 了 土地 整治 对 丘陵 田 块 利 用 产 
生 的 影响 。 孙 瑞 等 “ 选取 斑 块 个 数 、 平 均 斑 块 面 
耻 、 边 界 密度 、 斑 块 面 积 变异 系数 等 有 关 耕 地 细碎 
化 程度 的 景观 格局 指数 ， 分 析 了 有 效 表 征 耕 地 细碎 
化 变化 的 景观 格局 指数 的 适用 性 。 裴 欢 等 ' 通过 
景观 格局 指数 构建 了 脆弱 度 模 型 ， 并 以 景观 脆弱 
度 、 边 界 破碎 度 和 景观 类 型 破碎 度 构 建 了 耕地 景观 
生态 安全 评价 模型 。 吴 金华 等 ”利用 RRM 模型 研 
究 土 地 整治 风险 源 、 生 境 和 风险 受 体 之 间 的 相互 作 
用 关系 ， 对 土地 整治 划分 风险 等 级 ， 定 量 表达 出 土 
地 整治 的 生态 风险 值 。 此 外 ， 还 有 学 者 利用 数理 统 
计 方 法 对 景观 格局 与 耕地 特征 进行 相关 性 分 析 。 韩 
逸 等 …“ 通过 对 景观 格局 指数 与 耕地 质量 进行 相关 
性 分 析 ， 利 用 耕地 景观 格局 指数 反映 耕地 质量 。 张 
限 等 基于 全 局 回归 分 析 、 相 关 性 分 析 和 地 理 加 
权 分 析 人 研究 耕地 破碎 化 成 因 ， 发 现 耕 地 破碎 化 与 平 
均 坡 度 的 负 相 关 性 最 高 。 

目前 ， 这 些 研究 为 耕地 景观 格局 生态 风险 评价 
和 基于 景观 格局 的 相关 性 分 析 提 供 了 思路 和 方法 ， 
但 鲜 有 学 者 探究 景观 格局 生态 风险 与 地 形 因素 之 间 
的 关系 。 农 田 宜 机 化 整治 将 耕地 陡坡 改造 为 缓坡 
后 ， 耕 地 的 坡度 会 受到 影响 ， 从 而 造成 景观 格局 生 
态 风险 在 不 同 坡度 特征 上 呈现 出 分 布 差异 的 问 
题 中 。 评 价 农田 宜 机 化 整治 前 后 耕地 生态 风险 变 


化 情况 ， 以 及 探究 生态 风险 与 耕地 坡度 之 间 的 关系 
具有 重要 意义 和 应 用 价值 。 本 研究 以 重庆 市 潼 南 区 
20 个 县 为 研究 区 域 ， 基 于 2010 年 和 2020 年 土地 利 
用 数据 ， 采 用 ArcGIS 10.8 和 Excel 软 件 计 算 景 观 格 
局 指数 ， 通 过 烂 权 法 确定 各 指数 权重 并 构建 生态 风 
险 评 价 模型 。 对 2010 年 和 2020 年 的 数字 高 程 模型 
(Digital Elevation Model, DEM) 数据 进行 坡度 计算 
后 划分 坡度 等 级 ， 并 分 析 耕 地 坡度 不 同等 级 面积 
比 情况 。 从 景观 生态 学 角度 揭示 景观 格局 生态 风险 
时 序 和 空间 变化 特征 ， 并 对 景观 格局 生态 风险 与 耕 
地 坡度 不 同等 级 面积 占 比 进行 相关 性 分 析 ， 探 究 生 
态 风险 与 坡度 之 间 的 关系 ， 以 期 为 潼 南 区 农田 宜 机 
化 整治 提供 数据 参考 。 


1 数据 与 方法 
1.1 研究 区 域 概况 


潼 南 区 位 于 重庆 市 西北 部 ， 四 川 盆 地 中 部 偏 
东 ， 涪 江 中 下 游 ， 川 渝 交 界 处 ， 介 于 105°31'44"E~ 
106°00'20"E, 29°47'33"N~30°26'28" N Z lil], HUSH 
属于 川中 丘陵 区 ， 地 势 东 北部 和 南部 高 ， 中 部 低 ， 
区 内 大 多 数 地 区 海拔 在 250~-450 m 之 间 ， 地 势 总 高 差 
为 371 m; 坡度 小 于 25。 的 耕地 面积 占 全 区 耕地 面积 的 
90% 左 右 。 漳 南 区 属于 嘉陵 江水 系 ， 党 、 琼 两 江 从 西 
北向 东南 并 列 穿 过 区 境 。 漳 南 区 境 区 呈 狭 长 形 ， 东 西 
宽 46.9km， 南 北 长 72.1 km， 总 面积 1 584 km2， 辖 区 
共 3 个 街道 ，20 个 镇 。2015 年 ， 潼 南 区 率先 在 全 市 
进行 农田 宜 机 化 整治 ， 选 择 集 中 连 片 、 道 路 通达 、 
坡度 小 于 25° 的 地 区 ， 进 行 地 块 互联 互通 、 消 除 作 
业 死 角 、 优 化 地 块 布局 等 作业 ""。 截 至 2022 年 ， 
漳 南 区 完成 农田 宜 机 化 整治 800 pm”， 和 覆盖 22 个 街 
镇 。 根 据 各 镇 街 耕 地 数量 以 及 农田 宜 机 化 整治 情 
况 ， 本 研究 选取 潼 南 区 除 桂林 街道 、 梓 漳 街 道 、 大 
佛 街道 外 的 20 个 镇 为 研究 区 域 ， 高 程 图 如 网 1 
所 示 。 


1.2 数据 来 源 与 预 处 理 


W Pa DC 2015 一 2022 年 间 一 直 在 进行 农田 宜 机 
化 整治 。 为 探究 济南 区 农田 宜 机 化 整治 带 来 的 耕地 
景观 格局 生态 风险 在 时 序 和 空间 上 的 变化 ， 以 2015 
年 为 时 间 节 点 ， 选 择 整 治之 前 的 5 年 2010 年 ) 和 
整治 之 后 的 5 年 (2020 年 ) 作为 研究 时 期 ， 数 据 处 
理 均 基于 ArcGIS 10.8 #il Excel 软件 完成 。 两 个 时 期 
的 土地 利用 数据 来 自 地 理 监测 云 平台 网 站 (http:// 
www.dsac.cn/)， 从 中 提取 出 耕地 分 布 情况 的 栅 格 数 
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注 : QD 除 桂林 街道 .梓潼 街道 .大 佛 街 道 ; @ 该 图 基于 自然 资源 部 标 
准 地 图 服务 网 站 下 载 的 审 图 号 为 GS(2019)1822 号 标准 地 图 制作 , 底 
图 无 修改 。 
图 1 重庆 市 潼 南 区 高 程 图 
Fig. 1 Elevation map of Tongnan district of Chongqing city 


据 ， 用 于 分 析 耕 地 景观 格局 指数 。2010 年 数字 高 程 
模型 (Digital Elevation Model, DEM) 的 空间 分 辩 
率 为 12.53 m， 来 自 ALOS PALSAR 数据 集 ; 2020 年 
DEM 的 空间 分 辩 率 为 30 m， 来 自 SRTM1 V3.0 数 
据 集 (https:Wearthdata.nasa.gov/) ， 用 于 计算 耕地 坡 
度 。 中 国 行政 区 划 矢 量 数 据 来 自 自 然 资 源 部 标准 地 
图 服务 网 站 (http:Wbzdt.ch.mnr.gov.cn/) ， 用 于 裁剪 
土地 利用 栅 格 数据 和 DEM 栅 格 数据 。 
1.3 研究 方法 
1.3.1 景观 格局 指数 的 选取 

景观 格局 指数 可 以 定量 反映 景观 的 结构 组 成 和 
空间 配置 特性 ， 表 征 耕 地 的 地 块 规则 程度 、 细 碎 化 
规模 、 空 间 分 散 度 等 特征 "“"。 综 合 考 虑 各 指数 的 
含义 并 根据 农田 宜 机 化 整治 实际 作业 情况 ， 参 考 现 
有 的 土地 整治 以 及 景观 格局 生态 风险 评价 常用 指 
数 22， 选取 可 以 反映 耕地 形状 、 规 模 和 空间 分 布 
特性 的 6 个 景观 格局 指数 。 

潼 南 区 地 钥 以 丘陵 地 钥 为 主 ， 是 骨 型 的 浅 丘 平 
坝 区 域 ， 耕 地 地 势 较为 平坦 ， 基 本 以 条 田 和 格 田 为 
主 ， 其 空间 分 布 情 况 会 影响 耕地 使 用 效率 ， 所 以 选 


取 斑 块 密度 (Patch Density, PD) 和 分 离 度 (Divi- 
sion, D) 来 表征 耕地 地 块 的 分 布 情况 。PD 以 耕地 
斑 块 的 相对 大 小 和 密集 程度 来 反映 耕地 景观 破碎 
度 。PD 武 大， 耕地 斑 块 越 小 ， 越 难以 耕作 ， 耕 地 
破碎 化 程度 越 高 ““。D 表征 耕地 斑 块 的 分 离 程度 ， 
其 值 越 大 ， 说 明 耕 地 景观 越 分 散 ， 分 布 越 复 杂 ”。 

漳 南 区 丘陵 地 貌 中 的 耕地 地 块 零 星 、 分 散 ， 导 
致 农业 机 械 进 出 地 块 困难 、 无 法 作业 ， 严 重 制约 其 
农机 化 发 展 ， 阻 碍 漳 南 区 耕地 规模 化 发 展 ， 所 以 选 
取 分 维 数 (Fractal Dimension, FD) 和 边界 密度 
(Edge Density, ED) 表征 耕地 的 规模 化 程度 。FD 用 
于 测定 景观 斑 块 的 复杂 程度 ， 其 值 越 大 ， 耕 地 斑 块 
形状 的 复杂 程度 越 大 。ED 用 于 描述 耕地 景观 被 边 
界 分 割 的 程度 ， 直 接 反映 耕地 景观 的 破碎 程度 ”。 
此 外 ， 潼 南 区 丘陵 地 貌 中 的 耕地 形状 较为 复 
杂 ， 导 致 坡地 无 便道 可 通达 ， 而 且 农 业 机 械 在 地 块 
中 行进 、 转 弯 困 难 ， 使 得 耕地 更 加 难以 管理 ， 所 以 
选择 平均 形状 指数 (Mean Shape Index, MSI) 和 平 
均 斑 块 面积 (Mean Patch Size, MPS) 来 表征 耕地 
的 形状 特征 。MSI 由 形状 指数 (Landscape Shape 
Index, LSI) 进行 计算 ， 用 以 宏观 描述 耕地 形状 。 
LSI 接 近 1， 表 征 耕 地 形状 越 接近 正方 形 ， 更 易于 
管理 ; 反之 ，LSI 数 值 越 远 离 1， 则 表征 耕地 形状 
越 不 规则 ， 越 难以 管理 '”。MPS 用 于 宏观 描述 耕 
地 斑 块 面积 大 小 ， 斑 块 越 小 ， 耕 地 易 耕 性 越 差 。 相 
关 景 观 格 局 指数 生态 学 意义 及 计算 方法 如 表 1 
所 示 。 

13.2 ”生态 风险 评价 指数 的 构建 

为 定量 表达 耕地 景观 格局 变化 与 景观 生态 风险 
ZINA, Vi WUE ASE) 为 理论 基础 ， 基 于 
上 述 选取 的 6 个 耕地 景观 格局 指数 ， 构 建 景 观 干 扰 
度 指 数 和 景观 脆弱 度 指 数 ”， 并 根据 景观 干扰 度 
睛 数 和 景观 脆弱 度 指数 构建 生态 风险 评价 指数 ， 将 
耕地 景观 形状 、 规 模 、 空 间 特 性 定量 转化 为 景观 生 
态 风 险 。 

1) 景观 干扰 度 指 数 (Landscape Disturbance 
Index, UV) 。 在 漳 南 区 农田 宜 机 化 整治 过 程 中 ， 耕 地 
会 由 不 连续 、 复 杂 的 格局 变 为 互联 互通 、 相 对 简单 
的 格局 。 例 如 ， 农 田 宜 机 化 整治 可 以 减缓 耕地 坡 
度 、 消 除 作业 死角 、 拉 长 农业 机 械 行 进 路 线 、 增 加 
地 块 规模 等 ， 这 些 改造 情况 直接 影响 耕地 的 景观 格 
局 ， 能 够 减少 耕地 地 块 之 间 的 边界 线 和 过 渡 区 域 ， 
整体 上 提升 农田 的 规模 性 和 连续 性 ， 使 景观 格局 更 
加 统一 和 有 序 ， 所 以 选取 对 耕地 景观 格局 有 较 大 影 
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表 1 景观 格局 指数 生态 学 意义 及 计算 方法 
Table 1 Ecological significance and calculation method of landscape pattern index 
景观 格局 指数 生态 学 意义 计算 方法 

斑 块 密度 (PD ) 描述 景观 斑 块 破碎 程度 PD=n/A (1) 

分 离 度 (D) 描述 景观 斑 块 个 体 空间 分 离 程度 p=} j2 2) 

分 维 数 (FD) 描述 景观 斑 块 形状 复杂 程度 FD= 8) (3) 

边界 密度 (ED) 描述 景观 斑 块 被 分 割 程度 ED=s/A (4) 

平均 形状 指数 (MSI) 描述 景观 斑 块 的 平均 复杂 程度 Ms (5) 

平均 斑 块 面积 (MPS) 描述 景观 斑 块 大 小 MPS=A/n (6) 


TE: 其 中 ，7z 为 评价 单元 耕地 斑 块 总 数目 ， 个 ; APE Roce 


8 数 ; 计算 为 公式 (7): Ls 


啊 的 斑 块 密度 、 分 离 度 、 分 维 数 这 3 个 景观 格局 指 
数 构建 景观 干扰 度 指数 '””， 如 公式 (8) 所 示 。 

U=aX PD+b X D+c X FD (8) 

其 中 ,UU 为 景观 干扰 度 指 数 ; PD 为 斑 块 密度 ， 
km’; DD 为 分 离 度 ; FD 为 分 维 数 ; ar b, ca 
为 斑 块 密度 、 分 离 度 、 分 维 数 的 权重 。 

借鉴 已 有 研究 成 果 并 结合 研究 区 农田 宜 机 化 整 
治 的 实际 情况 ， 在 计算 干扰 度 指数 时 ， 认 为 斑 块 密 
度 对 景观 格局 的 影响 程度 最 大 ， 其 次 是 分 离 度 和 分 
维 数 ， 所 以 权重 a、b、c 分 别 赋 值 0.5、0.3、0.2， 
且 a 十 b 十 c=1 '*, 

2) 景观 脆弱 度 指数 (Landscape Vulnerability 
Index, V)。 景 观 脆弱 度 指数 表征 了 耕地 结构 的 脆弱 
性 ， 指 数 越 高 ， 生 态 系 统 越 不 稳定 。 根 据 各 指数 生 
态 学 含义 ,选取 与 耕地 脆弱 度 密切 相关 的 边缘 密 
度 、 平 均 形状 指数 、 平 均 斑 块 面积 这 3 个 景观 格局 
指数 构建 景观 脆弱 度 指数 ， 如 公式 (9) 所 示 : 

V=W ED + W,MSI-+W,MPS (9) 

其 中 , V 为 景观 脆弱 度 指数 ; ED 为 边缘 密度 ， 
km; MSI 为 平均 形状 指数 ;MPS 为 平均 斑 块 面积 ， 
km 个 ; Wi、W,、W; 分 别 为 边缘 密度 、 平 均 形状 
指数 、 平 均 斑 块 面积 的 权重 。 

3) 景观 脆弱 度 指数 权重 计算 。 综 合 现 有 景观 
格局 生态 风险 评价 中 的 权重 确定 方法 " "小 ， 相 比 于 
层次 分 析 法 等 主观 法 确定 权重 ， 烂 权 法 是 一 种 客观 
确定 权重 的 方法 ， 所 以 采用 该 方法 计算 边缘 密度 、 
平均 形状 指数 、 平 均 斑 块 面 积 这 3 个 景观 格局 指数 
的 权重 ， 从 而 提高 评价 结果 的 科学 性 。 

QD 归 一 化 计算 公式 。 由 于 景观 格局 指数 单位 不 
同 ， 存 在 正 向 指标 和 负 向 指标 的 问题 ， 在 计算 权重 


积 ，km ; s 为 评价 单元 斑 块 个 体 的 总 周 长 ，km; ZSZ 为 单个 斑 块 的 形状 
， 式 中 : /为 单个 斑 块 的 周 长 ，km; a 为 单个 斑 块 的 面积 ，km”。 


前 采用 极 差 标准 化 法 ”对 各 指数 进行 归 一 化 处 理 ， 
使 各 数据 具有 可 比 性 。 根 据 景观 格局 指数 对 耕地 景 
观 格局 的 影响 ， 将 边缘 密度 和 平均 形状 指数 划分 为 
负 向 指标 ， 平 均 斑 块 面积 划分 为 正 向 指标 。 负 向 指 
标 与 景观 脆弱 度 指 数 呈 负 相 关 ， 正 向 指标 与 景观 脆 
弱 度 指数 呈正 相关 。 归 一 化 计算 见 公 式 〈10) MA 
BC 11 Ye 

正 向 指标 : U= (X-X nn) / (Xina-X inn) (10) 

负 向 指标 : U= (Kimo X) / Xma Xion) (11) 

其 中 ， 忆 为 归 一 化 后 的 值 ; X; 为 第 i 个 县 耕地 
景观 格局 指数 的 原始 值 ，X;,,, 为 第 i 个 县 景观 格局 
指数 的 最 小 值 ; Xj 第 i 个 县 景观 格局 指数 的 最 
大 值 。 

OREM. EIIE TEA RE Ce) 计算 
方法 为 公式 (12). 

e=- [ln (n) ] > Un(D) (12) 

其 中 ,nn 为 评价 单元 个 数 , 个 。 

权重 (W,) 计算 方法 为 公式 (13): 

We (1-e,) X" e; (13) 

其 中 ，m 为 景观 格局 指数 个 数 ， 个 。 

计算 得 到 ED、MSI、MPS 的 权重 分 别 为 W= 
0.31、W,=0.45、W,=0.24。 

4) 生态 风险 评价 指数 (Ecological Risk Index, 
ERI)。 根 据 上 述 景观 干扰 度 指 数 和 景观 脆弱 度 指数 
构建 生态 风险 评价 指数 ， 表 征 评价 单元 生态 风险 程 
度 ， 指 数值 越 大 ， 表 明生 态 风险 越 高 ， 计 算 方 法 为 
公式 (14). 

ERI=UXV (14) 

其 中 ，ERI 为 生态 风险 评价 指数 ; 0 为 景观 干 
扰 度 指数 ; V 为 景观 脆弱 度 指数 。 
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13.3 ”生态 风险 值 与 坡度 相关 性 分 析 

根据 全 国 第 三 次 土地 调查 中 的 耕地 坡度 分 级 标 
准 ' 引 ,农田 宜 机 化 整治 尽量 选择 坡度 在 25° 以 下 区 
域 ， 其 中 15" 以 下 区 域 应 优先 考虑 进行 整治 ， 结 合 
研究 区 实际 耕地 坡度 地 形 特征 : 0°~2°* 和 >25° 坡 度 
耕地 面积 偏 少 ， 绝 大 部 分 耕地 坡度 在 2"~25"。 为 
探 明 农田 宜 机 化 整治 应 在 什么 坡度 范围 展开 工作 更 
加 适宜 ， 所 以 将 2°~25? 区 域 细 分 4 个 等 级 ,结合 研 
究 区 实际 耕地 坡度 地 形 特 征 ， 将 研究 区 耕地 坡度 按 
FA 0°~2°, 2°~5°, 5°~8°, 8°~ 15°, 15°~25°, >25° 
分 为 6 个 等 级 ， 分 别 记 为 一 级 、 二 级 、 三 级 、 四 
级 、 五 级 和 六 级 ， 并 分 别 计算 各 研究 区 域 不 同 坡度 
等 级 面积 占 比 。 基 于 SPSS 26 软 件 ， 对 坡度 等 级 面 
积 占 比 与 耕地 景观 格局 生态 风险 值 进行 相关 性 分 
析 ， 从 而 揭示 坡度 与 景观 生态 风险 之 间 的 关系 ， 进 
而 判断 农田 宜 机 化 整治 的 方向 。 


2 结果 与 分 析 


2.1 耕地 景观 格局 变化 分 析 
2.1.1 耕地 景观 生态 风险 时 序 特征 

根据 景观 格局 指数 计算 公式 ,计算 2010 年 和 
2020 年 各 评价 单元 的 耕地 斑 块 密度 、 分 离 度 、 分 维 
数 、 景 观 干扰 度 指数 、 边 界 密 度 、 平 均 形 状 指数 、 
平均 斑 块 面积 、 景 观 脆弱 度 指数 、 生 态 风 险 评价 指 
数 及 各 指数 的 平均 值 ， 结 果 如 表 2 所 示 。 从 2010 年 
到 2020 年 ， 斑 块 密度 、 分 离 度 、 分 维 数 平均 值 分 
别 由 0.75、1.76、0.35 降低 为 0.71、1.72、0.34， 表 
征 耕 地 的 破碎 度 降 低 ， 地 块 连 片 度 增加 。 干 扰 度 指 
数 平 均值 由 0.97 下 降 为 0.94， 耕 地 景观 格局 整体 抗 
干扰 能 力 增加 。 边 界 密度 平均 值 由 2.74 降低 为 
2.56， 平 均 斑 块 面积 平均 值 由 4.40 增 加 为 5.55， 耕 
地 被 分 割 程度 降低 ; 平均 形状 指数 平均 值 由 1.36 增 
加 为 1.53， 耕 地 斑 块 的 形状 更 加 不 规则 。 脆 弱 度 指 
数 平均 值 由 2.96 增 加 为 3.20， 耕 地 景观 格局 更 加 脆 
弱 。 生 态 风 险 评价 指数 由 3.10 下 降 为 3.01， 赫 地 景 
观 格局 生态 风险 降低 。 需 要 注意 的 是 ， 花 岩 镇 的 平 
均 形 状 指数 由 2010 年 的 5.10 增 加 为 2020 年 的 
27.25。 经 查询 原始 数据 得 知 ， 花 岩 镇 的 耕地 斑 块 
减少 ， 地 块 的 连通 性 增加 ,但 是 同时 增加 的 耕地 形 
状 的 不 规则 性 ， 所 以 会 出 现 数值 差异 较 大 的 现象 。 
漳 南 区 自 2015 年 开始 一 直 在 进行 农田 宜 机 化 整治 ， 
从 宏观 角度 看 ， 斑 块 密度 、 分 离 度 、 分 维 数 和 生态 
风险 的 平均 值 降低 ， 说 明 国家 实施 农田 宜 机 化 整治 
政策 有 显著 成 效 ， 降 低 了 耕地 破碎 度 ， 提 高 了 地 块 


之 间 的 连通 性 ， 并 提高 了 生态 安全 性 。 

表 2 2010 年 和 2020 年 重庆 市 江南 区 各 评价 单元 耕地 景观 
格局 指数 计算 结果 
Table 2 Results of cultivated land landscape pattern index cal- 
culation for each evaluation unit in Tongnan District of Chongq- 
ing city in 2010 and 2020 

评价 单元 年 份 PP D FD U ED MSI MPS V ERI 
© 2010 0.15 1.61 0.19 0.59 1.58 1.39 6.89 3.19 1.90 
2020 0.11 1.57 0.17 0.56 1.40 1.42 8.76 3.68 2.07 
2010 0.10 1.59 0.15 0.55 1.85 1.61 10.52 4.42 2.45 
2020 0.20 1.57 0.22 0.62 1.78 1.35 4.95 2.72 1.68 
_ 2010 0.67 1.86 0.41 0.98 3.13 1.28 1.48 2.18 2.14 
2020 0.48 1.82 0.35 0.86 2.89 1.31 2.09 2.26 1.93 
~ 2010 0.13 1.66 0.18 0.60 1.90 1.52 7.63 3.57 2.14 
` 2020 0.12 1.58 0.18 0.57 1.57 1.45 8.10 3.58 2.05 
~ 2010 0.11 1.64 0.16 0.58 1.71 1.42 9.38 3.99 2.31 
2020 0.13 1.61 0.18 0.59 1.57 1.32 7.47 3.36 1.97 
2010 0.20 1.51 0.22 0.59 1.80 1.36 5.10 3.49 2.07 
2020 0.04 1.42 0.10 0.46 1.54 4.60 27.25 9.47 4.40 
2010 0.15 1.65 0.19 0.61 1.87 1.40 6.59 3.23 1.97 
2020 0.09 1.57 0.15 0.55 1.53 1.45 10.54 4.29 2.35 
2010 0.17 1.64 0.21 0.62 1.95 1.43 5.76 2.99 1.85 
2020 0.33 1.71 0.29 0.73 2.23 1.35 3.04 2.29 1.68 
2010 0.21 1.51 0.23 0.60 1.55 1.31 4.82 2.52 1.52 
2020 0.14 1.49 0.19 0.55 1.46 1.34 7.29 3.23 1.78 
2010 0.33 1.84 0.29 0.78 2.90 1.40 3.00 2.57 2.00 
2020 0.24 1.78 0.24 0.70 2.52 1.39 4.22 2.78 1.95 
2010 0.24 1.61 0.25 0.65 2.21 1.42 4.16 2.79 1.82 
2020 0.18 1.59 0.21 0.61 2.13 2.12 5.51 3.01 1.84 
2010 0.19 1.67 0.22 0.64 1.83 1.30 5.22 2.75 1.76 
2020 0.14 1.59 0.19 0.58 1.48 1.29 7.19 3.21 1.87 
2010 0.14 1.72 0.18 0.62 2.26 1.54 7.33 3.61 2.24 
2020 0.32 1.68 0.28 0.72 2.06 1.34 3.09 2.25 1.63 
se SHU 2010 0.43 1.71 0.33 0.79 2.05 1.24 2.32 1.94 1.54 
2020 0.53 1.69 0.37 0.85 1.97 1.23 1.88 1.77 1.50 
2010 0.50 1.68 0.35 0.83 2.09 1.25 2.01 1.87 1.54 
2020 0.29 1.61 0.27 0.68 1.76 1.28 3.51 2.18 1.48 
2010 1.04 1.99 0.51 1.22 3.94 1.24 0.96 2.37 2.90 
2020 1.01 1.97 0.50 1.20 3.73 1.26 0.99 2.28 2.73 
2010 6.43 2.58 1.27 4.24 9.98 1.26 0.16 4.84 20.54 
2020 5.98 2.48 1.22 3.98 8.89 1.25 0.17 4.35 17.31 
2010 0.28 1.69 0.27 0.70 2.01 1.28 3.56 2.33 1.63 
2020 0.27 1.76 0.26 0.71 2.32 1.32 3,75 2.51 1.79 
2010 2.04 2.18 0.71 1.82 5.64 1.25 0.49 2.99 5.44 
2020 2.39 2.20 0.77 2.01 5.85 1.25 0.42 3.06 6.16 
2010 1.52 1.75 0.62 1.41 2.47 1.22 0.66 1.61 2.26 
2020 1.28 1.76 0.57 1.28 2.48 1.22 0.78 1.65 2.12 
2010 0.75 1.76 0.35 0.97 2.74 1.36 4.40 2.96 3.10 
2020 0.71 1.72 0.34 0.94 2.56 1.53 5.55 3.20 3.01 


2.1.2 ”生态 风险 空间 分 布 与 转化 特性 
根据 生态 风险 评价 指数 计算 得 到 的 耕地 景观 格 
局 生态 风险 值 ， 采 用 ArcGIS 10.8 软 件 对 各 评价 单 
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元 进行 赋值 ， 为 明显 区 分 各 生态 风险 等 级 的 不 同 ， 
根据 主观 分 级 法 "划分 为 5 个 风险 等 级 低 风险 


风险 (2.50<ERT) ， 得 到 2010 年 和 2020 年 的 耕地 景 
观 格 局 生态 风险 空间 分 布 图 (图 2) ， 并 计算 2010 


CERI<1.75) ， 较 低 风 险 (1.75<ERI<2.00), ， 中 风险 ”年 和 2020 年 两 个 时 期 各 生态 风险 等 级 的 面积 、 面 
(2.00<ERI<2.25), ， 较 高 风险 (2.25<ERI<2.50), ， 高 。 积 占 比 和 面积 变化 率 ( 表 3)。 
105°40'0"F 105°S0'0"F. 106°0'0"1 105°40'0"E 105°50'0"1 106°0'0"E 
= 一 = r 
é E z A Z 
Z E Z Z 
Z E Z Z 
z Zz 4 
105°40'0"E 105°50'0"E 106°0'0"E 105°40'0"T 105°S0'0'T 106°0'0"E 
a. 2010 b. 2020 
注 : 该 图 由 自然 资源 部 标准 地 图 服务 网 站 下 载 的 审 图 号 为 GS(2019)1822 号 标准 地 图 制作 , 底 图 无 修改 。 
图 2 2010 年 和 2020 年 重庆 市 潼 南 区 耕地 生态 风险 分 布 图 
Fig. 2 Ecological risk distribution maps of cultivated land in Tongnan District of Chongqing city in 2010 and 2020 
表 3 重庆 市 潼 南 区 生态 风险 面积 面积 占 比 、 变 化 率 
Table 3 Ecological risk area, area proportion, rate of change in Tongnan district of Chongqing city 
年 份 低 风 险 较 低 风 险 中 风险 较 高 风险 高 风险 
面积 /km” ”比重 /% ”面积 km” ”比重 /% 面积 km” 比重 /% ”面积 /km” 比重 /% 面积 km” ”比重 /% 
2010 310.68 22.32 482.16 34.62 217.77 15.64 214.66 15.42 167.19 12.01 
2020 330.79 23.75 511.90 36.76 271.08 19.47 84.36 6.06 193.43 13.96 
风险 面积 变化 大 小 /km 20.11 29.74 53.31 -130.30 26.24 
风险 面积 变化 率 /% 6.44 6.17 24.48 -60.70 16.30 


由 图 2 和 表 3 可 知 ，2010 年 和 2020 年 两 个 时 期 
高 风险 区 面积 比重 较 小 ， 面 积 增加 26.24 km ， 变 化 
率 为 16.30% ， 其 中 新 胜 镇 、 卧 佛 镇 、 五 桂 镇 风险 等 
级 未 变 ; 较 高 风险 区 面积 减少 130.30 km ， 变 化 率 
为 -60.70%， 其 中 玉溪 镇 、 古 溪 镇 、 宝 龙 镇 风险 等 
级 分 别 降低 为 中 风险 、 较 低 风 险 、 低 风险 ; 中 风险 
区 面积 增加 53.31 km*， 变 化 率 为 24.48%， 其 中 尝 
金 镇 风险 等 级 未 变 ， 太 安 镇 、 别 口 镇 风险 等 级 分 别 
降低 为 低 风险 、 较 低 风 险 ， 花 岩 镇 风险 等 级 转化 为 
高 风险 ; 较 低 风 险 区 面积 占 比 最 大 ，2010 年 与 
2020 年 较 低 风险 区 面积 比重 分 别 达 到 34.62% 与 
36.76%， 面 积 增加 29.74 km*， 变 化 率 为 6.17%， 其 
中 寿 桥 镇 、 双 江 镇 、 上 和 镇 风险 等 级 未 变 ， 米 心 镇 
风险 等 级 降低 为 低 风险 ， 柏 梓 镇 、 龙 形 镇 风险 等 级 


分 别 转化 为 中 风险 、 较 高 风险 ; 低 风险 区 面积 占 比 
仅 次 于 较 低 风险 区 ， 哪 种 面积 增加 20.11 km ， 变 化 
率 为 6.44%， 其 中 塘坝 镇 、 田 家 镇 风险 等 级 未 变 ， 
小 渡 镇 、 群 力 镇 风险 等 级 转化 为 较 低 风险 。 在 整个 
研究 期 间 ， 不 同等 级 的 风险 面积 变化 有 着 显著 性 的 
差异 ， 低 风险 、 较 低 风 险 、 中 风险 、 高 风险 区 域 范 
围 扩 张 ， 较 高 风险 区 域 范 围 缩小 ， 说明 风险 较 小 区 
域 在 扩张 ， 耕 地 整体 布局 上 风险 较 小 区 域 逐 渐 占 优 
势 ,农田 宜 机 化 整治 增加 了 地 块 作业 的 连通 性 ， 在 
小 地 块 合并 为 大 地 块 方面 取得 了 成 效 。 但 是 新 增 了 
高 风险 区 域 花 宕 镇 ， 经 查询 花 宕 镇 的 统计 年 鉴 ， 这 
可 能 与 其 近年 来 城镇 化 、 工 业 化 速度 加 快 有 关 。 
低 、 较 低 风 险 区 域 向 北部 蔓延 ， 中 风险 区 域 向 西部 
蔓延 ， 较 高 风险 区 域 仍 分 布 在 北部 ， 新 增 部 分 向 东 
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SEAS SUE, Fey AU Se DX STS EE) TE HBT 
ETB ay Ie] PY HE, FE TL, DT HS He Ae 
田 宜 机 化 整治 效果 较 好 ， 有 利于 开发 农业 资源 ; 南 
部 地 区 农田 宜 机 化 整治 效果 欠 佳 ， 后 续 农田 宜 机 化 
整治 工作 需 继续 优化 南部 地 区 耕地 景观 格局 ， 如 增 
加 地 块 的 连通 性 ， 消 除 地 块 间 的 过 渡 边 界 ， 以 保证 


70 
5 60 
与 s0 
is 和 
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20 
10 
0 
Zeek RRR EK REE RK ER HE HE KE BH 
# RO RM = RÆK ZSrReek ss 号 
E HE% EE] =Zo4oR RK FE E 
评价 单元 
s0 .2 wSRo mg%15° Ø15°-25° m25* 以 上 
a. 20104 


农业 的 可 持续 发 展 。 
2.2 生态 风险 与 坡度 等 级 分 析 
2.2.1 耕地 坡度 等 级 面积 占 比 分 析 

基于 Excel 软件 计算 各 评价 单元 耕地 不 同 坡度 
等 级 的 面积 占 比 ， 并 分 别 绘制 2010 年 和 2020 年 坡 
度 面积 占 比 柱状 图 ， 如 图 3 所 示 。 


WR Mm 


;FRERE SF BEE 


评价 单元 
m0~2° 25° m5%-8° wl5° W1525° 25° 以 上 


b. 2020 年 


图 3 2010 年 和 2020 年 重庆 市 潼 南 区 坡度 面积 占 比 柱状 图 
Fig. 3 Bar chart of slope area proportion in Tongnan district of Chongqing city in 2010 and 2020 


由 图 3 可 知 ，2010 年 和 2020 年 两 个 时 期 的 同一 
评价 单元 坡度 等 级 面积 占 比 变化 较 小 ， 但 不 同 坡度 
等 级 的 面积 占 比 有 明显 区 别 。 两 个 时 期 主要 以 四 级 
坡度 为 主 ， 面 积 占 比 介 于 23%~43% 之 间 ; 二 级 、 
三 级 坡度 面积 占 比 介 于 7%~35% 之 间 ， 且 同一 评 
价 单元 内 二 级 、 三 级 坡度 面积 占 比 大 小 接近 ， 两 者 
具有 相似 的 坡度 变化 趋势 ， 五 级 坡度 面积 占 比 介 于 
2%~38% 之 间 ， 不 同 评价 单元 面积 占 比 差异 最 大 ; 
一 级 、 二 级 坡度 面积 占 比 介 于 0.03%~13% 之 间 ， 
不 同 评价 单元 面积 占 比 差异 较 小 。 由 此 表明 ， 研 究 
区 耕地 坡度 大 多 处 于 2*~25"， 且 不 同 评价 单元 于 
此 坡度 区 间 的 地 形 起 伏 变化 较 大 ; 耕地 坡度 小 于 2° 
或 大 于 25° 区 域 面积 占 比 较 少 ， 且 不 同 评价 单元 于 
此 坡度 区 间 的 地 形 起 伏 变 化 较 小 。 

2.2.2 ”生态 风险 与 耕地 坡度 相关 性 
以 同一 坡度 等 级 下 的 各 评价 单元 坡度 面积 占 比 


表 4 重庆 市 潼 南 区 耕地 生态 风险 值 与 坡度 面积 占 比 的 相关 
Table 4 Correlation coefficient between ecological risk value of 


cultivated land and proportion of slope area in Tongnan district 


of Chongqing city 
坡度 等 级 相关 系数 

2010 年 2020 年 
一 级 -0.374 -0.471° 
二 级 -0.438° -0.534" 
三 级 -0.592" -0.677 © 
四 级 -0.320 -0.049 
五 级 0.609™ 0.729™ 
六 级 0.849™ 0.839™ 


注 : * 和 ** 分 别 表 示 在 P<0.05 水 平 显 著 相 关 和 P<0.01 水平 极 显著 
相关 。 

(0.592) > 二 级 (0.438) > 一 级 (0.374) > 四 级 
(0.320) ， 其 中 三 级 、 五 级 、 六 级 坡度 面积 占 比 与 
生态 风险 值 呈 极 显著 相关 ; 二 级 坡度 面积 占 比 与 生 


为 自 变量 ， 对 应 评价 单元 的 生态 风险 值 为 因 变量 ， 
基于 SPSS26 软件 对 不 同 坡度 等 级 面积 占 比 与 生态 
风险 值 进行 相关 性 分 析 ， 相 关系 数 结果 如 表 4 
所 示 。 

由 表 4 可 知 ， 不 同 坡度 等 级 面积 占 比 与 生态 风 
险 值 之 间 存 在 不 同 程度 的 相关 性 ，2010 年 和 2020 
年 两 个 时 期 的 一 级 、 二 级 、 三 级 、 四 级 坡度 面积 
比 与 生态 风险 值 呈 负 相 关 ， 五 级 、 六 级 坡度 面积 
比 与 生态 风险 值 呈正 相关 。2010 年 的 相关 系数 大 小 
依次 为 六 级 (0.849) > 五 级 (0.609) > 三 级 


态 风 险 值 呈 显 著 相关 ; 一 级 、 四 级 坡度 面积 占 比 与 
生态 风险 值 相关 性 不 显著 。 因 此 ， 在 0°~15° 坡 度 
范围 内 ， 生 态 风险 值 随 坡 度 面 积 占 比 的 增加 而 降 
低 ， 在 大 于 15? 坡 度 范 围 内 ， 生 态 风 险 值 随 坡度 面 
积 占 比 的 增加 而 增加 。2020 年 的 相关 系数 大 小 依次 
为 六 级 (0.839) > 五 级 (0.729) > 三 级 (0.677) > 
二 级 (0.534) > 一 级 (0.471) > 四 级 (0.049)， 其 
中 二 级 、 三 级 、 五 级 、 六 级 坡度 面积 占 比 与 生态 风 
险 值 呈 极 显著 相关 ; 一 级 坡度 面积 占 比 与 生态 风险 


Vol. 6, No. 3 


张 星 山 等 : 基于 景观 格局 的 


丘陵 山区 耕地 生态 风险 评价 65 


值 呈 显著 相关 ; 四 级 坡度 面积 占 比 与 生态 风险 值 相 
关 性 不 显著 。 在 0"~8" 坡 度 范围 内 ， 生 态 风险 值 随 
坡度 面积 占 比 的 增加 而 降低 ; 在 大 于 15° 坡 度 范 围 
内 ， 生 态 风险 值 随 坡度 面积 占 比 的 增加 而 增加 。 整 
体 来 看 ， 相 关系 数 随 坡度 的 增加 而 增加 ， 其 中 一 
级 、 二 级 、 三 级 、 五 级 、 六 级 坡度 面积 占 比 的 相关 
系数 呈现 出 严格 递增 趋势 ; 四 级 坡度 面积 占 比 的 相 
关系 数 最 小 。 两 个 时 期 六 级 坡度 面积 占 比 的 相关 系 
数 均 大 于 0.8， 说 明 耕 地 坡度 大 于 25? 时 对 生态 风险 
值 的 影响 最 大 ， 远 远大 于 耕地 坡度 处 于 8°~15° 范 
围 内 对 生态 风险 值 的 影响 。 此 外 ，2010 年 和 2020 
年 两 个 时 期 三 级 、 五 级 坡度 面积 占 比 的 相关 系数 有 
所 增加 ， 说 明 此 坡度 区 间 与 风险 值 的 相关 性 增强 。 
通过 以 上 生态 风险 与 耕地 坡度 相关 性 分 析 结 果 ， 在 
0? 一 15" 坡 度 范围 内 ， 为 平衡 低 坡度 的 生态 风险 ,在 
进行 农田 宜 机 化 整治 的 同时 ， 应 种 植 适宜 当地 条 件 
的 农作物 ， 以 提高 产量 和 质量 ; 在 大 于 15° 坡 度 范 
围 内 ， 为 平衡 高 坡度 的 生态 风险 ， 可 以 因地制宜 地 
确定 适宜 的 灌溉 方式 和 水 肥 一 体 化 管理 策略 ， 以 避 
免 土壤 侵蚀 和 水 资源 浪费 。 
3 讨论 与 结论 
3.1 讨论 

农田 宜 机 化 整治 主要 是 基于 耕地 景观 格局 进 
行 ， 耕 地 破碎 度 、 分 离 度 等 景观 格局 的 变化 会 影响 
耕地 的 质量 及 利用 效率 ， 所 以 采用 景观 格局 指数 构 
建 耕地 生态 风险 模型 可 以 评估 耕地 的 生态 风险 值 。 
此 外 ， 农 田 宜 机 化 整治 会 影响 整治 区 域 的 微 地 貌 、 
坡度 、 高 程 等 地 形 条 件 ， 而 坡度 大 小 为 影响 土地 侵 
蚀 的 主要 因素 ， 所 以 通过 分 析 耕 地 坡度 面积 占 比 与 
生态 风险 值 的 相关 性 可 以 确定 重点 整治 的 坡度 范 
围 ， 以 高 效 实施 农田 宜 机 化 整治 工作 。 

本 研究 通过 构建 生态 风险 指数 和 分 析 生 态 风险 
与 坡度 等 级 的 相关 性 ， 从 景观 格局 与 地 形 的 角度 对 
潼 南 区 耕地 生态 风险 进行 了 评价 ， 研 究 结果 如 下 。 

1) 根据 耕地 景观 格局 指数 计算 结果 ，2020 年 
相 较 于 2010 年 ， 研 究 区 域 斑 块 密度 、 分 离 度 、 分 
维 数 、 边 界 密度 平均 值 降低 ， 平 均 斑 块 面积 平均 值 
增加 ， 表 明 农 田 宜 机 化 整治 后 耕地 破碎 度 和 被 分 割 
的 程度 降低 ， 零 散 的 耕地 被 归并 后 规模 扩大 ， 耕 地 
连 片 度 增加 ， 更 加 利于 农业 机 械 集中 管理 和 规模 经 
营 ; 平均 形状 指数 增加 ， 耕 地 形状 趋 于 复杂 ， 农 业 
机 械 可 能 需要 考虑 进出 地 块 时 的 便利 性 问题 。 干扰 
度 指数 由 2010 年 的 0.97 下 降 到 2020 年 的 0.94， 耕 


地 整体 抗 干扰 能 力 增强 ; 脆弱 度 指 数 由 2.96 增 加 到 
3.20, 耕地 结构 更 加 脆弱 ， 这 可 能 与 近年 来 城镇 化 
速度 加 快 有 关 ， 部 分 耕地 的 原始 景观 结构 被 破坏 。 
生态 风险 评价 指数 由 3.10 下 降 到 3.01， 耕 地 生态 风 
险 值 降低 ， 农 田 宜 机 化 整治 有 效 改 善 了 耕地 景观 格 
局 ,一 定 程度 上 弥补 了 人 类 活动 对 耕地 景观 格局 的 
扰动 ， 提 升 了 农田 生态 系统 的 安全 性 。 

2) 耕地 景观 生态 风险 空间 分 布 及 转化 特性 显 
示 ，2010 年 和 2020 年 两 个 时 期 研究 区 主要 以 低 风 
险 和 较 低 风 险 区 域 为 主 。 低 风险 区 域 面积 增加 
6.44% ， 主 要 回 北 部 蔓延 ; 较 低 风险 区 域 面积 增加 
6.17% ， 主 要 向 中 东部 和 东南 部 草 延 ; 中 风险 区 域 
面积 增加 24.48 多 ， 主 要 向 西部 及 西北 部 草 延 ; 较 高 
风险 区 域 面积 减少 60.70%， 新 增 部 分 向 东北 部 草 
延 ; 高 风险 区 域 面积 增加 16.30%， 新 增 部 分 向 西北 
蔓延 。 整 体 来 看 ， 研 究 区 耕地 景观 生态 风险 趋 于 下 
降 ， 安 全 区 域 在 增加 。 

3) 由 耕地 生态 风险 与 坡度 相关 性 分 析 ， 
2010 年 和 2020 年 两 个 时 期 研究 区 耕地 坡度 主要 
以 2"~-25?" 为 主 ， 且 各 评价 单元 于 此 坡度 区 间 的 地 
形 起 伏 变 化 较 大 。 耕 地 坡度 小 于 1$*" 时 面积 占 比 与 
生态 风险 值 呈 负 相关 ， 耕 地 坡度 大 于 15° 时 面积 占 
比 与 生态 风险 值 呈正 相关 。2010 年 耕地 坡度 处 于 5° 
一 8*、15°~25°*、25° 以 上 时 与 生态 风险 值 呈 极 显著 
相关 ， 相 关系 数 分 别 为 0.592、0.609、0.849; 2020 
年 耕地 坡度 处 于 2°~5°、5°~8? 、15°~25° 、25° 以 
上 时 与 生态 风险 值 旦 极 显 著 相关 ， 相 关系 数 分 别 为 
0.534, 0.667, 0.729, 0.839, FAP HEHE UKE S°~8° 
和 15?"~25?" 的 相关 系数 有 所 增加 ， 坡 度 面 积 占 比 与 
生态 风险 的 相关 性 增强 。 根 据 相 关 资 料 ， 耕 地 坡度 
从 5$? 开 始 出 现 中 度 面 状 侵蚀 ，10?~15?" 出 现 强度 侵 
蚀 ， 从 降低 生态 风险 值 和 实际 地 形 情 况 考 虑 ， 农 田 
宜 机 化 整治 对 耕地 坡度 进行 整治 时 ， 应 重点 关注 漳 
南 区 南部 地 区 ， 并 主要 集中 于 5°~8°* 和 15°%~25° 区 
域 。 此 外 ,耕地 坡度 大 于 25°? 时 被 认为 不 适合 进行 
农业 生产 , 但 是 由 于 耕地 坡度 大 于 25° 时 与 生态 风 
险 高 度 相关 ， 可 以 适当 在 大 于 25° 的 小 范围 内 进行 
坡度 整治 以 减 小 坡度 。 

3.2 结论 


本 研究 通过 构建 生态 风险 评价 模型 ,揭示 了 重 
庆 市 漳 南 区 农田 宜 机 化 整治 前 后 生态 风险 时 序 空间 
变化 特征 ， 并 基于 数理 统计 原理 ， 对 生态 风险 与 坡 
度 进行 相关 性 分 析 ， 探 究 了 生态 风险 与 坡度 的 关 


66 智慧 


业 ( 中 英文 ) Smart Agriculture 


Vol. 6, No.3 


系 。 通 过 比较 2010 年 和 2020 年 两 个 时 期 的 生态 风 
险 值 ， 农 田 宜 机 化 整治 降低 了 耕地 生态 风险 ， 能 提 
高 耕地 生态 风险 的 安全 性 ; 耕地 坡度 面积 占 比 与 生 
态 风险 值 具有 相关 性 ， 从 潼 南 区 实际 地 形 和 降低 生 
态 风 险 角度 来 讲 ， 后 续 农 田 宜 机 化 整治 工作 应 重点 
关注 漳 南 区 南部 耕地 处 于 $?"~8" 和 15?"~25" 区 域 。 

本 研究 未 考虑 耕地 内 部 因素 ， 如 土壤 有 机 质 含 
量 、 土 壤 板结 程度 等 对 生态 风险 的 影响 ， 今 后 可 以 
通过 加 入 土壤 质量 相关 因素 更 加 全 面 分 析 农 田 宜 机 
化 整治 对 耕地 的 生态 风险 的 改变 情况 。 


利益 冲突 声明 : 本 研究 不 存在 研究 者 以 及 与 公开 
研究 成 果 有 关 的 利益 冲突 。 
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Abstract: 

[Objective] Farmland consolidation for agricultural mechanization in hilly and mountainous areas can alter the landscape pattern, ele- 
vation, slope and microgeomorphology of cultivated land. It is of great significance to assess the ecological risk of cultivated land to 
provide data reference for the subsequent farmland consolidation for agricultural mechanization. This study aims to assess the ecologi- 
cal risk of cultivated land before and after farmland consolidation for agricultural mechanization in hilly and mountainous areas, and 
to explore the relationship between cultivated land ecological risk and cultivated land slope. 

[Methods] Twenty counties in Tongnan district of Chongqing city was selected as the assessment units. Based on the land use data in 
2010 and 2020 as two periods, ArcGIS 10.8 and Excel software were used to calculate landscape pattern indices. The weights for each 
index were determined by entropy weight method, and an ecological risk assessment model was constructed, which was used to reveal 
the temporal and spatial change characteristics of ecological risk. Based on the principle of mathematical statistics, the correlation 
analysis between cultivated land ecological risk and cultivated land slope was carried out, which aimed to explore the relationship be- 
tween cultivated land ecological risk and cultivated land slope. 

[Results and Discussions] Comparing to 2010, patch density (PD), division (D), fractal dimension (FD), and edge density (ED) of 
cultivated land all decreased in 2020, while meant Patch Size (MPS) increased, indicating an increase in the contiguity of cultivated 
land. The mean shape index (MSI) of cultivated land increased, indicating that the shape of cultivated land tended to be complicated. 
The landscape disturbance index (U) decreased from 0.97 to 0.94, indicating that the overall resistance to disturbances in cultivated 
land has increased. The landscape vulnerability index (V) increased from 2.96 to 3.20, indicating that the structure of cultivated land 
become more fragile. The ecological risk value of cultivated land decreased from 3.10 to 3.01, indicating the farmland consolidation 
for agricultural mechanization effectively improved the landscape pattern of cultivated land and enhanced the safety of the agricultural 
ecosystem. During the two periods, the ecological risk areas were primarily composed of low-risk and relatively low-risk zones. The 
area of low-risk zones increased by 6.44%, mainly expanding towards the northern part, while the area of relatively low-risk zones in- 
creased by 6.17%, primarily spreading towards the central-eastern and southeastern part. The area of moderate-risk zones increased by 
24.4%, mainly extending towards the western and northwestern part, while the area of relatively high-risk zones decreased by 60.70%, 
with some new additions spreading towards the northeastern part. The area of high-risk zones increased by 16.30%, with some new ad- 
ditions extending towards the northwest part. Overall, the ecological safety zones of cultivated relatively increased. The cultivated 
land slope was primarily concentrated in the range of 2° to 25°. On the one hand, when the cultivated land slope was less than 15°, the 
proportion of the slope area was negatively correlated with the ecological risk value. On the other hand, when the slope was above 
15°, the proportion of the slope area was positively correlated with the ecological risk value. In 2010, there was a highly significant 
correlation between the proportion of slope area and ecological risk value for cultivated land slope within the ranges of 5° to 8°, 15° to 
25°, and above 25°, with corresponding correlation coefficients of 0.592, 0.609, and 0.849, respectively. In 2020, there was a highly 
significant correlation between the proportion of slope area and ecological risk value for cultivated land slope within the ranges of 2° 
to 5°, 5° to 8°, 15° to 25°, and above 25°, with corresponding correlation coefficients of 0.534, 0.667, 0.729, and 0.839, respectively. 
[Conclusions] The assessment of cultivated land ecological risk in Tongnan district of Chongqing city before and after the farmland 
consolidation for agricultural mechanization, as well as the analysis of the correlation between ecological risk and cultivated land 
slope, demonstrate that the farmland consolidation for agricultural mechanization can reduce cultivated land ecological risk, and the 
proportion of cultivated land slope can be an important basis for precision guidance in the farmland consolidation for agricultural 
mechanization. Considering the occurrence of moderate sheet erosion from a slope of 5° and intense erosion from a slope of 10° to 
15°, and taking into account the reduction of ecological risk value and the actual topographic conditions, the subsequent farmland con- 
solidation for agricultural mechanization in Tongnan district should focus on areas with cultivated land slope ranging from 5° to 8° 
and 15° to 25°. 
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